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O ■ Abstract 
(N ■ 

' We discuss inhomogeneous cosmological models which satisfy the Copernican principle. We con- 

, struct some inhomogeneous cosmological models starting from the ansatz that the all the observers in the 

' models view an isotropic cosmic microwave background. We discuss multi-fluid models, and illustrate 

O . how more general inhomogeneous models may be derived, both in General Relativity and in scalar- 

tensor theories of gravity. Thus we illustrate that the cosmological principle, the assumption that the 

5-H . Universe we live in is spatially homogeneous, does not necessarily follow from the Copernican principle 

and the high isotropy of the cosmic microwave background. We also present some new conformally flat 

, two-fluid solutions of Einstein's field equations. 

•1— I . 

X: 

■ 1 Introduction 

The standard model of cosmology rests on several fundamental assumptions. As with any theoretical model 
of a physical system, it is crucial that these assumptions are identified and tested wherever possible, in order 
that the proposed model be considered acceptable. The standard model of cosmology is founded on the 
(perturbed) spatially homogeneous and isotropic cosmological models of Friedman-Lemaitre-Robertson- 
Walker (FLRW), which are derived from the cosmological principle. The cosmological principle may be 
taken to state that the universe is spatially homogeneous This is a strong assumption; considerably 
stronger than the Copernican principle which says that we are not at a special location in the Universe. 
Regardless of which principle one cares to take, when these are combined with assumed perfect isotropy 
about ourselves, on all scales, we arrive at the FLRW models. If we assume such perfect isotropy about us 
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without a 'mediocrity' principle, then we must be at a center of symmetry. Obviously, the properties of a 
spatially homogeneous universe can be radically different from a spherically symmetric one, for example, 
and it is therefore of fundamental importance to find some method to test the assumption of homogeneity, 
and to identify what exactly happens if any of the assumptions are relaxed. This is particularly important, 
since studying inhomogeneous models will allow us to identify some possible tests of non-homogeneity in 
the universe. 

What evidence do we have that the universe is so isotropic about us? Obviously the strongest and most 
important piece of evidence for this is the extremely high isotropy of the cosmic microwave background 
(CMB) which is isotropic to one part in 10^. The question is: In what context can we infer spatial homo- 
geneity from our observations of the CMB? Without the Copemican principle (or something similar) the 
answer is obviously not - we may be located at the 'center' of the universe and thus see the CMB isotrop- 
ically distributed about us, whereas all other observers would not have such a unique and privileged view. 
However, if we assume the Copernican principle (i.e., we assume that all or 'most' observers in the universe 
see the CMB to be as isotropic as we see it) can we infer homogeneity on the basis of the CMB alone? 

The first attempt to answer this question resulted in a theorem by Ehlers, Geren and Sachs [2| (hereafter, 
EGS) which states that if all observers in an expanding dust universe see an isotropic radiation field then 
that spacetime is homogeneous and isotropic (and therefore FLRW). This can trivially be generalised to the 
case of a geodesic and barotropic perfect fluid Q. The 'isotropic radiation field' is implicitly identified with 
the CMB. However, as has been emphasised recently [T 5 1, the resulting spacetime will be FLRW only if 
the matter content is of perfect fluid form, and the observers geodesic and irrotational. This work has been 
extended [4| to include inhomogeneous universe models with non-geodesic observers. That is, inhomoge- 
neous spacetimes have been found which also allow every observer to see an isotropic CMB. It has also been 
shown that a significant subset of these models are consistent with other observational constraints, regardless 
of observer position |6|. This means that these models are consistent with observations on all scales even 
when the Copernican principle is taken into account - and yet the models are significantly inhomogeneous. 
However, a problem of these models is that the non-geodesic motion of the observers prohibits a barotropic 
equation of state for the matter (although the models admit a thermodynamic scheme). There have been 
recent developments along this line, where in 1 7 1 a realistic multi-fluid form of the matter was proposed 
(although the models used were slightly different from those in |4| which allow an isotropic radiation field). 
Alternatively, a fractal distribution may be more appropriate ||8l. Other recent work concerning rotating and 
anisotropic cosmological models also supports these results ||9l ^1 ; other non-FLRW models exist 

which admit an isotropic radiation field. 

The purpose of this paper is to discuss more general cosmological models which allow an isotropic 
radiation field. Specifically, we wish to discuss models with 'realistic' matter in order to demonstrate that 
there exist physically viable inhomogeneous cosmological models which will allow an isotropic radiation 
field by construction but are not FLRW. The recent supemovae data imply an accelerated expansion rate in 
the universe: within the standard model this implies some sort of negative pressure, be it a cosmological 
constant or quintessence or some other type of exotic matter. Therefore, we will consider here not just 
traditional barotropic perfect fluid matter, but more exotic forms, such as scalar fields and varying A models. 

Spacetimes which allow all observers the view of an isotropic CMB must satisfy the 'isotropic radiation 
field theorem'. The isotropic radiation field theorem may be derived from the Einstein-Boltzmann equations 
for photons in a curved spacetime. It is easy to show from the multipole expansions of |3| that a spacetime 
with an isotropic radiation field must have the velocity field, of the photons being shearfree and obeying 

Ua = ^aQ, e = 3Q, (1) 
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where Q is a function of the energy density of the radiation field. Any observers travehng on this congruence 
will observe the isotropic radiation. This velocity field is also a conformal Killing vector of the spacetime. 
In fact, a spacetime admitting an isotropic radiation field must be conformally stationary, and we use this 
fact to construct some 'generalised-EGS' spacetimes. 

In the following we show how some irrotational multifiuid spacetimes may be constructed, which sat- 
isfy Q. To this end, we consider two non-comoving perfect fluids, which can be interacting or non- 
interacting, and may or may not admit barotropic equations of state for the fluids. The fluids are chosen 
to be non-comoving to allow for energy flux and anisotropic pressures in the energy momentum tensor; 
otherwise the models may be written as a single perfect fluid, which are the models studied in |4||6||. In a 
similar vein we consider models with a perfect fluid and scalar field; in contrast to usual work on mixtures 
of this kind, we allow for the case where the scalar field has a spatial gradient relative to the perfect fluid, 
which we take to be comoving with the isotropic radiation (test) field. In both cases we take one of the 
perfect fluids to be 'comoving' with the radiation field; that is one of the fluid velocities will be (parallel to) 
the timelike conformal Killing vector of the spacetime. We then consider some models with non-zero heat 
flux, but zero anisotropic pressure, which were previously considered in lil3J . which may be interpreted as 
'quintessence' models with varying A and energy flux. 

The case of non-zero rotation has also been considered elsewhere, and simple expanding and rotating 
spacetimes with plausible matter in which the observers could measure an isotropic CMB have been con- 
structed fTOlfmiQlim. This may be considered as a counter-example to numerous claims that the rotation 
of the universe may be constrained by observations of the CMB alone; such results make additional assump- 
tions of the matter present and the velocity field we follow in the universe. 

The upshot of all this is to emphasise that the high isotropy of the CMB when combined with the Coper- 
nican principle is simply not enough to draw conclusions about the spatial homogeneity of our universe. 

2 Spacetimes admitting an isotropic radiation field 

We are interested in spacetimes in which the high isotropy of the CMB is permissible for every observer. 
In the particular case where we have a model in which all observers on some congruence u"" see an exactly 
isotropic radiation field, then this velocity field has two important properties: 



isotropic radiation field are conformally stationary, with the velocity fields of the two (conformally related) 
spacetimes parallel - see the appendix. Now, if we were to assume that these observers measured only dust, 
then that spacetime must be FLRW - the origional EGS theorem |2|. 

In this paper, for simplicity, we also restrict our attention to the irrotational case. This may be justified by 
the following considerations. If part of the matter consists of a conserved comoving barotropic perfect fluid 
other than radiation, or for geodesic motion with any matter source, it follows from (Q that the expansion or 
the rotation must be zero. For a conserved barotropic perfect fluid, we have Uq = Va^, and p'9 = cj), where 
(j) = — J dp/ {fJ.{p) +p), and p' = dp/djj.; so, r]abc^'^V'^{Q — 0) = 2(| —p')Ouja = 0. For geodesic motion, 
Vahc^^^'^Q = \G^a = 0. However, rotating universes which allow an isotropic radiation field have been 
found and discussed in some detail - see fTO*. '303121. 

In this case the metric can then take the form 




(2) 



Writing iia - ^9ua = 



V aQ we see that the first condition is equivalent to Q. Spacetimes admitting an 




(3) 
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where Hap{x'^') can be diagonalised. \f Q = Q{t) then the acceleration is zero and we recover the models 
studied by Coley and MacManusj 15|; indeed, even in this case (i.e., the acceleration-free case) it follows 
that there are physically viable spacetimes that are not FLRW. 

In order to find irrotational spacetimes with an isotropic CMB, we can simply compute the Einstein 
tensor of and equate with the matter we desire. In the appendix, we discuss this computation further. 

2.1 The energy momentum tensor of multiple fluids 

In general, the energy momentum tensor. Tab for any spacetime may be decomposed with respect to the 
velocity field in the following covariant manner: 

Tab = P-UaUb + phab + 2q(^aUb) + T^ab- (4) 

This decomposition allows us to make the physical interpretations that p = u°-v}'Tab is the energy density, 
p = ^h°'^Tab the isotropic pressure, qa = —haU^Tbc the energy or heat flux and itab = T(^ab) the anisotropic 
pressure or stress. All these quantities are interpreted by an observer traveling on the congruence. 

2.1.1 perfect fluids 

Consider the energy-momentum tensor due to two non-comoving perfect fluids; 

Tab = fJ-lUaUb + Plhab + fJ'2UaUb + P2hab, (5) 

where fii are the energy densities of the fluids in each comoving frame, and the pi's are their respective 
pressures. The velocity field of the second congruence may be written as a Lorentz boost of the first; 

r = 7(n'^ + t^"), 7= \ , v^ua = 0. (6) 
If we write this as one fluid with respect to the u'^ congruence, then Tab has the form of @, with components 



p 


= W^u'Tab 




p 


1 i,abrp 
— 3"- -Lab 


= PI+V2 + (/"2 + P2 


qa = 


-hS'Tbc 


= 7^{f^2+P2)Va, 



T^ab = Ti^ab) = l^{f^2+P2)V(^aVb)- (7) 

Thus we see that the first fluid will experience an energy flux due to the second fluid passing through their 
frame (provided /X2 + P2 7^ 0). 

2.1.2 perfect fluid plus scalar field 

In general the energy-momentum tensor of a scalar field may be written 

= 4>,a4>,b - gab (^(^.c^'^ + V{(l>)) , (8) 

which, when a velocity field is specified, takes the form 

Ti = i^^UaUb + Va4>^b4> - '^4'U(a^b)<^ " 9ab (|Ve</>V^(/. - + V {c\>)) . (9) 
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Thus, if we add to this scalar field a perfect fluid (perfect with respect to this congruence), then we find 
that the total or mean matter variables become 



p = p + P0 = p + i(/>'-iv,(/.v^<^-y(0), 

= qt = -4>^a(t>, 

TTafe = T^tb = ^{a<P^b)4>- (10) 

Note that if Va4> = then formally the total fluid takes the form of a single perfect fluid. 

We may demand that satisfies the Klein-Gordon equation, which may be derived from the energy 
conservation equation for the scalar field, V"T^^ = 0, 

= V,V> = n'^V,0 + V'^V.^ -4>-9^. (1 1) 

ocp 

However, if there is an interaction between the scalar field and some other matter, then this equation may 
not hold; for example, we could have a scalar field decaying into physical matter in which case the energy 
of the scalar field will not be conserved independently. 

2.2 Scalar Tensor Theories of Gravity 

Scalar-tensor theories, in which a long-range scalar field combined with a tensor field mediate the gravita- 
tional interaction, are standard alternatives to general relativity. The original motivation for these theories 
was to incorporate a varying gravitational constant into GR to account for alleged discrepancies between 
observations and weak-field predictions of GR. A special case of the scalar tensor theories, known as the 
Brans-Dicke theory of gravity (BDT 1 16 1) (with a constant uq parameter), was the original of these theories. 
Scalar-tensor theories occur as the low-energy limit in supergravity theories from string theory 1 17 1 and 
other higher-dimensional gravity theories [18|. 

Recently the recovery of the EGS theorem in scalar tensor theories was given fi9ll ; geodesic observers 
in a scalar tensor theory of gravity observing isotropic radiation must be in a FLRW universe. We mention 
here that inhomogeneous spacetimes are possible however if the geodesic assumption is dropped. The field 
equations, obtained by varying the BD action with respect to the metric and the field (/>, are 

Gab = ^Tab + ^ (Va4>Vb(l) " ^5afe V^-Z-V^^) + (/."^(VaVb^ - (7abVcV"0) - i<?abC/(0). (12) 

(3 + 2a;)VaV> = SttT - VawVV + ^ (13) 

where the energy-momentum tensor of the matter, T"*, may take any of the usual desired forms. In the 
scalar-tensor gravity theories the principle of equivalence is guaranteed by requiring that all matter fields are 
minimally coupled to the metric gab- Thus energy-momentum is conserved: 

V"T„b = 0. (14) 

It is known that scalar-tensor theories can be rewritten in the conformally related 'Einstein' frame |20|, so 
that the models are formally equivalent to GR coupled to a scalar field. Therefore, scalar-tensor theories 
may be incorporated here as a special case of the scalar field in GR - see Sec. l2. 1.2l 
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3 Some solutions 



Rather than provide an exhaustive study of multi-fluid solutions, we will present some example of how 
such solutions may be derived. This is in keeping with our aim of illustrating the existence of 'realistic' 
inhomogeneous cosmological solutions which satisfy the Copemican principle. 



3.1 Two perfect fluids 



For simplicity, we restrict ourselves to the case where Q{t,x°') = Q{t,r) in comoving polar coordinates 
and we also restrict ourselves to the spherically symmetric case, Hrr = Hqq = H^^ = exp 25 (r, 9,<j)). In 
this case the energy flux relative to the comoving u"" frame is 



qr = 2e ^ {Q^tr - Q,tQ,r) = 7^(/"2 +P2)Vr, 



^ Vg = = 0. 



(15) 



Because v"- has one non-vanishing component, this implies from Q that iTab must be diagonal. Calculating 
the components of iiab, we find that B can only be a function of r alone: 



TTrr = - -{B" + 2Q") + {B' + 2Q'f 

2 qr ' 

We may also calculate the mean energy density and pressure; 



2Q'^ + -{B' + 2Q'] 
r 



Hence we find 



(16) 



(17) 



-2{Q+B) 



2{B" + Q") + {B' + Q'f + -{B' + Q' 



Ml + Ai2 + l^V^ {[12 + P2) 



,-2(Q+B) 



2(B" + Q") + + 5Q'2 ^ ^ ^ 2Q') 

r 



Pi + P2 + hl^v^ {fJ'2 + P2) ■ 



(18) 



As yet we have only determined v (and 7), and we have three equations relating four functions, Hi,pi 
(Eqs. (fTSl . and the remaining freedom from vf^r and Qr). In principle we have the freedom to specify one 
more equation relating the four free functions. The most obvious restrictions are barotropic equations of state 
for the two fluids, pi = Pi{ni), or separate energy conservation for the two fluids, V"T^j = 0. However, 
if Q and B are specified then only one of these types of conditions may be used in general. We can use 
the freedom in the two metric functions to allow us to use both conditions if we choose. For simphcity, we 
consider the case of the two fluids obeying linear equations of state, pi = Wifj,i. 



3.1.1 Example: pi = Wi^ii, with B = 

We have four free functions, Q in, and Vr together with four equations; specifying two equations of state 
is then sufficient to close the system. We may solve Eqs. (I15t and (IT6l for /i2 and v,,- as functions of Q 
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(and derivatives of). Substituting these into (fTSb we get two equations for ^i; requiring equality leads to a 
horrendous equation for Q; 



= {{2r'^ + 2wir'^)Qrr'^ + {{—wir"^ — r"^ + ^W2Wir'^ + 2>W2r'^)Qr'^ + {&W2Wir + 2wir 
+ 4 W2r)Qr + (-2 W2r'^ -2r^)Qu + {-W2r^ - 3 W2Wir^ - - 3 wir'^)Qt'^)Qrr 
+ (— 3Ti;2'^^ — Wir'^ — (iw2Wir'^ — r'^)Qr^ + {—7w2r — 5wir — 3r — 9w2Wir)Qr^ 
+ ((2 tuar^ + 2 r2)Qii + (3 W2r'^ + + tf^ir^ + 3 W2Wir'^)Qt'^ - 4:W2 - 4:Wi - 2 - 6 W2Wi)Qr'^ 
+ ((2 W2r + 2 r)Q44 + (3 wir + ^i;2^' + 3 W2Wir + r)(5t^ + ( 

+ 2w2r^Qtr^ -2wir'^Qtr^)/rwi{-Qr-rQr^ + rQrr)il + W2). (19) 
A solution of the form Q = alnt + blnr exists, provided we choose 

_2 l + W2^ + 2w2 ^20) 

—8 -0)2101 — 2W2 — <iWi + 3 W2^ — 2wi^ + 2) W2^Wi — 6 W2Wi^ — 1 ' 

h = -2 3^2u;i + 2u;2 + u;i ^^l) 
1 + 3 W2 + ^^1 + 3 W2W1 

The only other physical constraint is that < 1, for all t and r. In the case where the first fluid is dust, 
wi = 0, this requires that — | < W2 < — g- This also ensures that the fluids become comoving at late times. 

Thus we have demonstrated the existence of inhomogeneous two barotropic fluid solutions of the field 
equations which allow the existence of isotropic radiation. There are clearly much more general solutions 
than we have presented here, our solution being a very special case. 

3.2 Perfect fluid plus scalar field revisited. 

As before, we restrict ourselves to the spherically symmetric case where Q{t,x") = Q{t,r) and Hj.r = 
Hq0 = Hfj,^ = exp(2i?(r, 9, (/))). In this case the energy flux relative to the comoving u"- frame is 

Qr = 2e-^{Q^tr - QtQr) = -<P^A, Qd = = ^ ^{t, r). (22) 

Because is only a function of time and the spatial coordinate r, this implies that Ttab must be diagonal. 
Calculating the components of itab, we find that B can only be a function of r alone. We find that 

= Ttrr = l[-iB" + 2Q") + 2Q'^ + B'^ + 4B'Q' + -{2Q' + B')]. (23) 
3 3 r 

Hence we find 

We may also calculate the mean energy density and pressure; 

/2 = ^e-^QQ'^-e-''^'^+''\2{B" +Q") + {B' + Q'f + -{B' + Q')] 

r 

= /i + ^cA' + ^V.cAV^ + (24) 
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P = -^[3e-2«[Q2 + 2Q] - e-2(Q+iJ)[4g" + 2B" + AB'Q' + + B'^ + -{B' + 2Q')]] 
6 r 

= p + ^<^2-^V.</.V'-</)-F(</.), (25) 

We may impose the restriction of a barotropic equation of state for the perfect fluid p{i^) and an energy 
conservation law for the fluid V^Tab = 0. We choose the equation of state p = w^. The energy conservation 
equations for the perfect fluid read 

^i,t + ^Q,t{l + w)^i = {) (26) 
w^i' + Q'{1 + w)^x = 0, (27) 
which when integrated imply that Q must have the form 

Q{t,r) = a{t)+P{r). (28) 

It follows that a specific form for the comoving energy will be 

^ = exp[-(l + ii;)(3a + /3/tt;)]. (29) 
To simphfy things we will assume that 5 = 0, so that the metric now becomes 

9a, = e^«W+^'^('-)r?„,; (30) 
which impUes that 6 = 3a,te '5, iir = P', and 

? r .. .o 9 .1 

(31) 

Qr = -2e-^atP'. (32) 

Hence we find that 
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2 ^ Qr _ 2 2^/2 1 



The scalar field wave equation is 



-4> + - 2Q<j> + (2Q' + B')cp'e-'^ + ^e'^^ = ^e^«. (34) 

r d(f> 

We will assume a solution of the form 

= $(t) + ^(r), (35) 
which allows the derivation of the following two equations from the scalar field equation when V = Q: 

4> + 2d$ = C (36) 

+ 2(/3' + -)*' = C, (37) 
r 

where C is a constant. We can rewrite the equations for heat conduction and anisotropic pressure as 

= 2e-^d/3' (38) 
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^'2 = _2/3" + 2/3'V-/3' (39) 
r 

Hence we compute the following differential equations for a and 13: 

a + a^ = Ae" (40) 

/?" + + 2^ = Be^ (41) 
r 

These equations have non-trivial solutions implying a non-FLRW cosmology (the FLRW limit is recovered 
when [3 = 0). Hence we have shown that spacetimes with a barotropic perfect fluid and a non-comoving 
scalar field exist which allow an isotropic radiation field for all observers, which are non-FLRW. Clearly 
there are a huge number of solutions meeting this criteria; we have demonstrated existence in this simplest of 
cases. These new solutions could play an important role in cosmology, for example as a new generalisation 
of quintessence. 



4 Inhomogeneous quintessential cosmologies 

Simple multifluid models can be constructed by introducing, together with a barotropic fluid, a varying A(t) 
term as the second fluid: 

/) = ^ + A, p = (7 - 1)^ - A (42) 

where p, p are "total" energy density and pressure obtained from the field equations. The varying A term can 
be interpreted as the asymptotic state of a scalar field associated with a quintessence dominated scenario, 
coexisting with a material fluid described by ^. We show in this section that such an interpretation is 
compatible with the asymptotic properties of a class of simple models that allow an isotropic radiation field. 

The simplest models satisfying the EGS criterion and compatible with the decomposition (l42l are char- 
acterized by the conformally FLRW metric (a particular case of ^) given by 

.2 -dt2 + dr2 + r2(d02^sin2 0d(^2) ,,^2 
ds^ = ^ 1-^, <^ = a{t) + b{ty (43) 

whose source is the imperfect fluid 

Tab = P UaUb +pKb + '^q{aUb) (44) 

where = ^6"^ t and Qa = qr a> while p and p must comply with Applying the field equations for 
(l43l and (ISJ leads to 

Stt 7/x = 2 (a,4t + 26 + 6 ^r^) (a + br^) (45) 

87rg, = -4b,tr, 8tt q = Stt Ig'^^aQbl^^^ = 4\bt\r{a + br^) (46) 
8^ 7A = L4 (t) r^ + 2L2{t)r^ + Lq (t) (47) 

where 

Uit) = 3blj-2b,ttb, L2{t) = 3jatb,t-2b^-aub-b,tta, Lo{t) = 87 (a^, + 4a6) - 2a (26 + a «) 

(48) 
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If we demand that A be only a time-dependent function, we obtain 

L2{t) = 0, Uit) = 0, Stt A = Lo{t)/j (49) 

which yields differential equations that determine a, b for a given 7(t), and the definition of A(t). 

Since we ai^e interested in an asymptotic regime that assumes a slowly varying 7(t), we shall consider a 
constant 7. The general solution of the system ( l49t in this case is 

6 = 6o^"^^ a = al^-2^ + a2^-^^'^-^^^(^-^)^ 7^2/3, v=l/{3j-2) (50) 

where 60 1 «!) '3^2 are arbitrary integration constants. Since we are assuming fi characterizes a material fluid 
(baryons plus photons and possibly CDM), we have that 1 < 7 < 2 and hence u in ( 1501 is always positive. 
Inserting (|50|) into (|43ll, ^5^, (|46l), and (|49ll leads to 



1 _ 1 _ t^"'" 



(51) 



_ 12 [(27 - 5/3)60^3 + (fli + bpr^) t + (37 - 1) 02] [(-bo/3) f + (ai + fepr^) t + 02] 

8.,r = lr Svr , = 8|^ol IW) +^(^1 + boQ t + a2kr ^^^^ 

^ ^ ^ (-8/3)6g t6 + 12aibo + l6a2K + Saj 

^4(37-1)1/ ^ 

Because of the apparent (and coordinate dependent) resemblance of (l43l to a spatially flat FLRW, it is 
tempting to assume that these two metrics have common geometric features. For example, it is evident from 
d52l . d53l and (l54t that fi, A and q diverge as t — > for 1/ > and 1 < 7 < 2, hence we can identify 
t = as the locus of a "big bang" singularity, analogous to the FLRW big bang. In spatially flat FLRW 
spacetimes it is always possible to assume that the coordinate range is given by < t < cxd and < r < cxd, 
so that t — > 00 and r ^ 00 mark asymptotic future infinities in the timelike, null and spacelike directions. 
However, for the models under consideration the coordinate domain is necessarily restricted by the extra 
condition that the conformal factor !/<!> be a bounded function. Also, proper time along the worldlines of 
comoving observers is r = J dt/^ evaluated for fixed {r,9,(j)), and so a sufficient condition for having 
T ^ CXD occurs if the conformal factor l/# diverges, even if it does so for finite values of the coordinates 
t,r. From dSTT i. this occurs for all ai, 02, 60,7, since the equation (— 6o/3)t^ + (ai + bor^)t + 02 = 
always has real roots in the coordinate domain t > 0, r > 0, defining the hypersurface 



B = [t,r{t),e,(l)], r{t) 



(l/3)5ot^ - ait - a2 
bot 



(55) 



which can be represented as a parametric curve in the t,r plane. If one or both of ai, 02 is zero, the 
boundary B persists, though its parametrization in the t, r plane is simpler than d55l . The only exception is 
if 60 = 0, whence the solutions trivially reduce to FLRW. Therefore, this feature is inherent to the models 
characterized by and d44l i. 

The fact that and Y = r<I>~^ diverge at (l55l means that B marks a spacetime boundary beyond 
which the spacetime manifold cannot be extended. An asymptotic past/future is then defined as the coor- 
dinate values marked by B which are reached by causal curves, either comoving observers (r = const.) or 
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radial null geodesies: (v = t + r, w = t — r). From dSTt . it is straightforward to prove that r ^ oo holds 
as B is reached by future and past directed worldlines of comoving obervers. Also the affine parameter of 
radial null geodesies diverge at spacetime points marked by B. The coordinate domain of definition is then 
restricted by |<1>| > and depends on the signs of the constants oi, 02, bo, specifying the form of B in the 
plane t, r. ^From the various numerical values for these constants, we eliminate all those cases in which the 
evolution of the comoving observers occurs between two branches of B. The remaining cases display the 
two types of evolution classified below: 

Case (i). If oi, 02, 60 negative, B lies in the infinite past of all observers evolving towards their infinite 
future as t — > 00, and we have a null infinity analogous to that of a FLRW cosmology (the infinite past 
is then marked by B). Using null coordinates the asymptotic limit along outgoing radial null geodesies 
w ^ 00 is given by 

g^^^l66|(37j^ ^2(3,-4).^ STrg- ^^2(3,-4).^ ^ ^ _^ ^2(3,-4). (56) 

3 3 

so that a regular null infinity requires 7 > 4/3 (otherwise, the affine parameter has a finite limit as t ^ 00 
and this locus marks a null singularity). For a heat conducting shear-free fluid the weak energy condition 
requires: p + p = 7// > 2(7 1 21 1, a relation that is satisfied bv the asvmptotic forms d56l onlv for 2 / 3 < 7 < 
1. Hence, this case is unphysical. 

Case (ii) If ai, 02, 60 are positive, then all worldlines of comoving observers start their evolution at t = 
(big bang) and evolve towards their infinite future at B. From d52t . (l53l and d54t . we have 

q 2u^\bQ\t'^r 



fi (27 - 5/3) bot^ + (ai + 6or2) t + (37 - 1) 02 

A _ 1^2 [(-(8/3) bot^ + 12 aif^ + 16 02*^) bo + 3a|] 

^ ~ 12 [(27 - 5/3) 60^^ + (ai + bor^) t + {3j - 1) 02] [{-bo/3) + (oi + bor^) t + 02] 

so that near t = we obtain 



q If^r A V 



2 



H (27-5/3)t' n 4(37-1) 

while at the boundary B we have 



(57) 



(58) 



(59) 



^^0, g^O, p^A, p^-A, A^A^=A(t^) (60) 
^v^^Mtl^B . _ Abotl[{ai-2borl)t^+2a2] 



B 2botl+a2 ' " tt(=^7-i)- 



B 

where , are related by The limits ( 15^ indicate (for 1 < 7 < 2) that the models are matter 

dominated at the big bang (0 < A ^ fi), evolving towards a A dominated future at B (represented by dMT l. 
Notice that the asymptotic future state at B can be de Sitter (or anti de Sitter) and is not an asymptotically 
homogeneous state since A^ depends on position (for each observer a different constant value). For = 
(in (I55I ) we have that A^ > 0. As grows along B, Ag < 0, which implies that the total energy density 
p = /i + A becomes negative asymptotically. Hence, for the physical reasons, we exclude coordinate values 
r > fg, where satisfies A^ (f^ ) = 0. 

The behavior of q is compatible with the energy conditions, since q ^ fir/t <^ holds all along the 
evolution, near the big bang and near B. It is still necessary to find an adequate physical interpretation for 
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this term, whether as a heat flux or as a kinetic term associated with a velocity field or the dipole of a kinetic 
theory distribution However, since we are interested mainly in the asymptotic stage near B, as long as 
r is sufficiently small we will have q <^ fJ. and could consider g as a residual term. 

Another feature of the models, absent in FLRW spacetimes, is the fact that fi and q diverge as r — > oo 
along hypersurfaces of constant t that do not intersect B, marking a point singularity at spacelike infinity. 
This feature is also present for perfect fluid sources of (l43l . see |23 1. Note that models similar to those ex- 
amined here were considered recently flT\ ; however, the existence of the boundary ( 1551 was not considered 
in the asymptotic study of those models. 

The models discussed in this section illustrate how even simple inhomogeneous spacetimes have a much 
richer geometrical structure that heavily constrains their physical applicability. Pending a reasonable physi- 
cal interpretation for qa and provided we exclude sufficiently lai^ge values of r, these solutions are inhomo- 
geneous models that comply with the EGS criterion and describe a A dominated scenario usually associated 
with the "quintessence" field L24..25..26J . 

5 Discussion 

In this paper we have proven the existence and examined the physical viability of a number of spacetimes 
which have been constructed to allow an isotropic radiation field. Since these inhomogeneous spacetimes 
satisfy the Copernican principle (as far as the CMB is concerned), the question of finding methods of testing 
the cosmological principle, and thus observationally testing whether the universe is in fact an FLRW model, 
arises. 

More precisely, it has been shown here and elsewhere |6l that inhomogeneous universe models with 
non-geodesic observers obey the EGS criterion. That is, inhomogeneous spacetimes have been found which 
allow every observer to see an isotropic CMB. It has also been shown that a significant subset of these 
models are consistent with other observational constraints, and hence these models are consistent with ob- 
servations even when the Copernican principle is taken into account - and yet the models are not spatially 
homogeneous |6|. A potential problem with these particular models is that the non-geodesic motion of the 
observers prohibits a barotropic equation of state for perfect fluid matter. However, we have shown here 
that more general and physically viable cosmological models (with realistic matter) allow an isotropic radi- 
ation field. In particular, irrotational multi-fluid spacetimes have been constructed which satisfy Q- These 
cosmologies include two non-comoving perfect fluids, which can be interacting or non-interacting, and may 
or may not admit barotropic equations of state for the fluids. The fluids are chosen to be non-comoving to 
allow for energy flux and anisotropic pressures in the energy momentum tensor (otherwise the models may 
be written as a single perfect fluid and correspond to the models studied in |4 6|). Even in the acceleration- 
free case there are examples of spacetimes that are not FLRW (TSi . Similarly, models with a perfect fluid 
and scalar field can be constructed in which the scalar field can have a spatial gradient relative to the perfect 
fluid, which is taken to be comoving with the isotropic radiation field. As a particular example, a class 
of shear-free spherically symmetric, inhomogeneous (quintessential) cosmologies whose source is a heat 
conducting fluid and a scalar field were considered in detail. 

One of our key assumptions has been zero rotation. It has been shown that rotating spacetimes which 
allow an isotropic radiation field may also be constructed |'l0l[nn3'9l. 

Other recent work also supports these conclusions. In the fundamental EGS theorem 13, and here, it is 
assumed that all fundamental observers measure the CMB temperature to be exactly isotropic during a time 
interval / (defined by tE < t < to, where is the time of last scattering and to is the time of observation). 
Under this assumption the theorem then asserts that the universe is exactly an FLRW model during this time 
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interval. However, the EGS theorem cannot be used to conclude that the physical universe is close to an 
FLRW model since the CMB temperature can only be observed at one instant of time on a cosmological 
scale. Hence it is of interest to ask what restrictions, if any, can be placed on the anisotropy in the rate of 
expansion, assuming that all fundamental observers measure the CMB temperature to be exactly isotropic at 
some instant of time to only. On the basis of continuity, it can then be argued that all fundamental observers 
will measure the CMB temperature to be almost isotropic in some time interval of time of length 6 centered 
on to- This time interval could, however, be much shorter than the time interval /. However, in |12| it 
was shown that for a given time to, there is a class of locally rotationally symmetric non-tilted dust Bianchi 
type VIII spatially homogeneous cosmological models such that at to the CMB temperature is measured 
to be isotropic by all fundamental observers, even though the overall expansion of the universe is highly 
anisotropic at to- 

In addition, the EGS theorem is of course not directly applicable to the real universe since the CMB 
temperature is not exactly isotropic. This result has consequently been generalized by |27| to the "almost 
EGS theorem", which states that if aU fundamental observers measure the CMB temperature to be almost 
isotropic during some time interval in an expanding universe, then the universe is described by an almost 
FLRW model during this time interval. The dimensionless shear parameter and the Weyl parameter were in- 
troduced in [28 1 . Since the Weyl curvature tensor is related to time derivatives of the shear tensor, restricting 
the shear parameter to be small does not guarantee that the Weyl parameter is small. Therefore a necessary 
condition for the universe to be close to an FLRW model is that both of the shear and Weyl parameters must 
be small. In the almost-EGS theorems the dimensionless time and spatial derivatives of the multipoles are 
assumed to bounded by the multipoles themselves |27|. If this assumption is not satisfied, then the CMB 
temperature observations do not impose upper bounds on the shear and Weyl parameters, and hence do not 
establish that the universe is close to FLRW. In ["29^1 a class of spatially homogeneous non-tilted Bianchi 
type VIIq dust models in which the CMB is treated as a test field or a non-interacting radiation fluid was 
studied. To obtain the present CMB temperature pattern, the photon energies were integrated numerically 
along the null geodesies that connect points of emission on the surface of last scattering with the event of 
observation at the present time. Wainwright et al. [30 1 then showed that the shear parameter tends to zero but 
the Weyl parameter does not tend to zero at late times in these models. In other words, although the models 
isotropize as regards the shear, the Weyl curvature remains dynamically significant. A variety of numerical 
simulations to calculate anisotropy patterns of the CMB temperature in Bianchi VIIq models were explicitly 
performed [29| to demonstrate that there exists cosmological models that are not close to any FLRW model 
even though the temperature of the CMB is almost isotropic in the sense that the observational bounds on 
the quadrupole and octupole are satisfied. 

It is clear, then, that the Copernican principle when combined with our observations of the CMB does 
not imply the cosmological principle: that the universe is homogeneous and isotropic. Recent work on this 
suggests that the assumptions in the EGS (dust observers) and almost-EGS (small gradients of CMB mul- 
tipoles and dust observers) theorems are crucial to the conclusions; weakening of any of these assumptions 
appears to negate the theorems almost entirely. It is therefore important to ask what observations we need 
to test the cosmological principle? 

Assuming for the moment that the assumptions of the almost EGS theorem actually hold in our universe, 
then one method by which the cosmological principle may be tested is as follows: if we can observe the CMB 
as seen by some other observers, then we can immediately confirm or reject the cosmological principle. That 
is, if we find the CMB is as isotropic around these other observers as we see it around us, we may conclude 
that we live in a homogeneous universe. There is, in fact, a physical method by which we can observe the 
CMB as seen by other observers. It consists of light from the CMB being scattered by hot gas in galaxy 
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clusters in such a way as to allow us to observe the anisotropy of the CMB as seen by that particular galaxy. 
This is known as the Sunyaev-Zel'dovich (SZ) effect OTl . and has been suggested by a number of authors 
as a possible means to test the Copernican and cosmological principles (3T\. 

This will only work of course if the initial assumptions of the EGS theorem apply in our actual Universe. 
The original EGS theorem relies on the observers in the universe being well described by a dust fluid - i.e., 
they are geodesic. Indeed, in the almost EGS theorem, it was necessary to demand dust observers to first 
order - i.e., more general matter was only allowed at second order. 

Any possible attempt to verify the cosmological principle by using methods such as the SZ effect above, 
or any other method which relies on making observations of the CMB from other locations, will fail. If 
we were living in a universe found here or in |4|, for example, we would see exactly the same effect: all 
observations of the CMB around other observers would be as isotropic as the standard homogeneous FLRW 
models. It follows, therefore, that the high isotropy of the CMB can never be used, on its own, to show our 
universe is nearly homogeneous. 

On the other hand, we can test the cosmological principle using the SZ effect and methods hke it pro- 
vided we can show definitively that our universe is made of a dust-like fluid and that we travel on geodesies. 
Recent observations suggesting quintessential matter making up a significant part of the energy density of 
the universe throws this standard assumption into question: there is no a priori reason why this matter - 
whatever it turns out to be - should be homogeneous (many dark matter theories allow equations of state 
which are not dust also). 

The geodesic assumption can be tested to some degree by local observations. If we look closely at the 
recession velocity of galaxies close to us, then we can detect a dipole moment in the relative velocities of 
these galaxies. This deviation from the hnear Hubble law is usually attributed to our local random motion 
with respect to the local expansion rate of the universe, caused by our gravitational infall into the Great 
Attractor. Acceleration will also leave its mark on the linear Hubble law also as a dipole distribution in 
the direction of the CMB dipole, but as a dipole which grows linearly with distance [4|. With sufficiently 
accurate knowledge of the local distribution of galaxies, the two effects can be disentangled as our bulk 
gravitational motion affects the Hubble law irrespective of distance. Current knowledge of our local group 
motion will only provide relatively weak constrains on acceleration (roughly u/Hq ~ 0.1; compare this 
with a /Hi ~ 10"^ from almost-EGS CMB observations [33 1). 

If we assume, for arguments sake, that the results of such a study reveal that the local dipole is entirely 
accounted for by our peculiar velocity, then what can we say about the spatial homogeneity of our universe? 
If the acceleration around our location is small, then we may assert the Copernican principle and ascertain 
that all observers in the universe follow geodesies; therefore we may apply the EGS theorem (assuming the 
SZ measurements measure small enough anisotropics of the CMB around other observers) and deduce that 
the cosmological principle is a valid assumption. This would yield tremendous support for our faith in the 
standard model. 

Of course, acceleration is just one possible inhomogeneity which causes problems: others are rotation, 
Weyl curvature, and anisotropy of the energy momentum tensor. These all leave their mark to varying 
degrees to the anisotropy of the magnitude-redshift and number-count-redshift relations, but at higher order 
in redshift (if expanded in a power series in redshift) than acceleration, so are much harder to detect. 

What should be clear from this argument is that the cosmological principle is simply untestable. Even if 
SZ measurements reveal the Copernican principle to be true, accurate determination of anisotropics of the 
magnitude-redshift relation at high redshift are out of the question for the foreseeable future. 
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A. Appendix. The energy momentum tensor of conformally related spacetimes 



A conformal transformation is an angle preserving transformation that changes lengths and volumes. The 
importance of these types of transformations lies in the fact that, under a conformal transformation, the 
causal structure of the spacetime is preserved. The Weyl tensor, C^^^^^, is invariant, so that a conformal 
transformation will introduce no tidal forces or gravitational waves; that is, a conformal transformation will 
only introduce 'non-gravitational' forces and matter into the new spacetime (by changing Rah and thus the 
matter tensor Tab via Einstein's equations). 

We will discuss conformal transformations and their 1+3 splitting here as it is a useful tool for con- 
structing new spacetimes from old, especially if one is after spacetimes with an isotropic radiation field. 
The conformal transformation is defined by 

(62) 



gah 



e^Ua 



where Q > is an arbitrary function, u"" is a velocity vector with respect to gab- gabU°'u^ = UaU' 



,a 



-1; and 
= -1. 



is the conformally related (parallel) velocity vector, and is normalised with respect to gab'- gabu"'u^ - 
The covariant derivative of any one-form field Va transforms as 

VaVfe = VaVb - 2Q(^aVb) + OabQ^Vc, (63) 

where Qa = Q,a = ^aQ — Qua- The expansion (9 = divn), acceleration (ii"' = u^Vbu"-), rotation 
{uja = — |curl Ua), and shear (Uab = ^(a'^J-b)) of the two velocity congruences are related by: 

e = e^{9-3Q) 



VaQ 



Ua 



Ua 

^ab = e'^aab- (64) 

The equation for the acceleration corrects equation (6.14) of l34l . These show that a conformal transforma- 
tion may induce acceleration and expansion into the new spacetime, but not shear or rotation: in particular, a 
conformally flat model must have vanishing shear and rotation fT W|. With respect to gah, a dot denotes time 
differentiation along the fluid flow. Fab = u^^cFab, and Va is the spatial derivative projected orthogonal to 
the flow lines, V cFab = h^c^'^a^^b^d.Fef, where hab = gab + UaUb is the usual projection tensor. 
The Einstein tensor transforms as 

->2 



Gab — Gab — aQb — '^QaQb + gab 2V cQ'^ — Q 



(65) 



where Gab is the Einstein tensor of gab, and 
time and space derivatives: 



QaQ'^- For clarity, we decompose derivatives of Q into 



VbQa 



Qa = ^aQ - QUa 
UaUb [Q - ifVcQ 



2u 



(a 



-Vb)Q + ^9Vb)Q + (o-fe)c + Vb)dc^' 



+ k^ab 



V^VcQ-Q9 -QcTab + VaVbQ. 



(66) 
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We also write Tab = Gab, and Tab = Gab as general fluids, both with respect to u"". 

Gab = P-Uaitb + PKb + "^QiaUb) + T^ab, (67) 
Gab = llUaUb + phab + 2qi^aUb) + T^ab] (68) 

where {fi, p, qa, VTab}, and {fi, p, Qa, TTab} are the energy density, isotropic pressure, heat flux, and anisotropic 
pressure of Gab and Gab respectively. 

We can decompose Gab given by (1551 into the fluid variables in (|5HJl by using ( I66t in the following 
covariant manner: 

H = u^vf'Gab = e^^fi -3q(^Q- '^9^ - 2VaV''Q + VaQV^Q, (69) 
P = Ih'^'Gab = e^'^p +(q-^9^Q-2Q + ^V^V'^Q - ^V„QV»Q + lu^V'^Q, 

Qa = -u''G(^a)b=e'^qa + 2Q(VaQ-Ua)+2h^(VbQy +4VabcUJ^^''Q (70) 
VTafe = Gl^ab) =T^ab + 2Qaab-2Vl^a^b)Q-'^^{aQ^b)Q- (71) 

In the energy flux equation we used the identity 

VaQ = (VbQ) ' -QUa + ^9VaQ + a.^VfoQ + r]abcOJ^V'Q- (72) 
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